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We study the elementary excitations of a model Hamiltonian for the 7r-electrons in poly-diacetylene 
chains. In these materials, the bare band gap is only half the size of the observed single-particle gap 
and the binding energy of the exciton of 0.5 eV amounts to 20% of the single-particle gap. Therefore, 
exchange and correlations due to the long-range Coulomb interaction require a numerically exact 
treatment which we carry out using the density-matrix renormalization group (DMRG) method. 
Employing both the Hubbard-Ohno potential and the screened potential in one dimension, we 
reproduce the experimental results for the binding energy of the singlet exciton and its polarizability. 
Our results indicate that there are optically dark states below the singlet exciton, in agreement with 
experiment. In addition, we find a weakly bound second exciton with a binding energy of 0.1 eV. The 
energies in the triplet sector do not match the experimental data quantitatively, probably because 
we do not include polaronic relaxation effects. 

PACS numbers: 71.20.Rv, 71.10.Fd, 78.30.Jw, 78.20.Bh 



I. INTRODUCTION 
A. Experimental observations 

1. Structure 

Poly-diacetylenes (PDAs) are prototypical quasi one- 
dimensional materials. ^'^ Their monomer building unit 
is comprised of four carbon atoms. The four outer elec- 
trons of each carbon atom are sp^ hybridized. Three of 
them form cr-bonds. The cr-bonds are between neighbor- 
ing carbon atoms on the chain and to covalent ligands 
R and R' , which are several A long and differ for var- 
ious members of the PDA family. The fourth electron 
is delocalized over the carbon backbone in a molecular 
TT-orbital. 



the single and double bonds alternately form angles of 
^p^ = 120° and Lp2 ^ 240° degrees. 

Very long undistorted polymer chains have been built 
starting from a monomer single crystal so that the chains 
are perfectly ordered,^ and single polymer chains diluted 
in their monomer matrix have even been prepared and 
studiedi^ In PDAs, exciton-polaritons have been gener- 
ated that have been shown to be coherent over tens of 
micrometers, i.e., several ten thousand monomer units.— 
Consequently, the opto-electronic properties of the PDAs 
result from the electrons' mutual interaction and their 
interaction with the periodic lattice potential, while the 
influence of disorder is negligible. 



2. Optical properties 



FIG. 1: Lewis structure of a poly-diacetylene unit cell. 



The resulting Lewis structure is shown in Fig. [TJ The 
four carbon atoms in the unit cell are linked by a triple 
bond, a single bond, a double bond, and a single bond. 
The atomic distances are = 1.4 A, = 1.3 A, and rt — 
1.2 A for the single (S), double (D), and triple bonds (T), 
respectively. The chain of atoms is not perfectly straight; 



PDAs are insulators; the gap for single-particle excita- 
tion is A ~ 2.4 eV. Band structure calculations estimate 
the bare band gap to be Abare ~ 1.2 eV, i.e., electronic 
exchange and correlations account for half of the single- 
particle gap. Moreover, the optical gap for the primary 
exciton is Aopt — 1.9 eV, so that the exciton binding en- 
ergy Abind — 0.5 eV is about 20% of the band gap. Due 
to the restricted geometry, the electron-electron interac- 
tion must be treated accurately for the calculation of the 
optical properties of the PDAs. 

Lattice effects complicate the analysis of the spectra 
of PDAs in two ways. First, the primary exciton with 
excitation energy Aopt is accompanied by phonon side- 
bands that result from oscillations of the double and 
triple bonds. These signals dominate the optical exci- 
tations below the band gap. Second, some PDA single 
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crystals such as TCDIJ^ exist in two different conforma- 
tions, which have exciton energies Abiuo = 2.0 eV ('blue 
chains') and Aiod — 2.4 eV ('red chains'), respectively. 
Therefore, it is not easy to disentangle the effects of the 
electron-electron interaction from those of the electron- 
lattice interaction; the latter gives rise to resonance shifts 
of several tenths of an electron volt. 



B. Theoretical approaches 

1. Extended Wannier theory 

In order to describe the optical excitations in PDAs, 
two approaches have been taken. The first approach 
starts from an ab-initio density-functional theory calcu- 
lation of the bare band structure in local-density ap- 
proximation (LDA), which is then supplemented by an 
approximate treatment of the residual electron-electron 
interaction, e.g., the GW approximation for the single- 
particle bands and the Bethe-Salpeter equation (BSE) 
for the excitons (LDA-|-GW^-|-BSE).^ Actual calcula- 
tions for the PDAs often omit the GW step (Wannier 
theory) 4 Within this approach, a number of experimen- 
tal data could be reproduced, e.g., the exciton binding 
energy and its polarizability. Unfortunately, the theory 
does not predict optically dark states below the exciton 
resonance. 



2. Model calculations 

The second approach to a theoretical description of 
the primary excitations in polymers starts from a many- 
particle model Hamiltonian that describes only the tt- 
electrons and their mutual interaction. The parameters 
for the kinetic energy of the electrons are taken from LDA 
calculations, and the Coulomb interaction is approxi- 
mated in various ways, e.g., with the Ohno parametriza- 
tion^ of the Pariser-Parr-Pople potential.^" With the 
help of the DMRG methodjii the ground state and ele- 
mentary excitations for such models can be calculated for 
large chains with a very high accuracy. In this way, the 
electron-electron interaction is treated without resorting 
to any approximations. Unfortunately, specific calcula- 
tions for the PDAai^ have only had limited success in 
meeting the experimental test. 



3. Simplifying assumptions and outline of this work 

In this work, we consider structures which start with 
a triple bond and end with a double bond, i.e., we con- 
sider the sequences [(TSDS)m_iTSD] of = 4to carbon 
atoms with m triple bonds and m double bonds. We 
perform a DMRG study of a many-body model Hamilto- 
nian for the TT-electrons which uses Barford and Bursill's 



parametrization for the band-structure partiS, but em- 
ploys the screened Coulomb potential in one dimension. 
The essential difference between the screened potential 
and the Ohno potential is that the local interaction part 
is larger in the screened potential. Because of this, the 
Hubbard-Ohno potential leads to a good description of 
the low-energy singlet excitations in PDAs. 

We do not include lattice relaxation for the single- 
particle and optical excitations. As has been shown 
by Barford and co-workers the relaxation energies for 
triplet excitons in trans-polyacetylene can be as large as 
0.4 eV. The parameter sets we will use below apply to 
the rigid-lattice situation and will change slightly when 
lattice relaxation is taken into account properly. 

The different varieties of PDA differ in their ligands, 
which introduce local potentials on the carbon atoms 
with double bonds. Therefore, the excitation energies 
of the PDAs differ by a few tenths of an electron volt. In 
this work, we ignore the ligand effects and only consider 
a prototypical case. 

II. MODEL FOR POLY-DIACETYLENE 

A. Operators in second quantization 

1. Kinetic energy operator, current operator, and dipole 
operator 

In this work we will restrict ourselves to the descrip- 
tion of the TT electrons because they dominate the optical 
response of the poly-diacetylenes for energies TujJ < 3 eV. 
The motion of the electrons is described by the operator 
for the kinetic energy, 

f^-^ti (c+^Q+i,^ + cf^^/ci^^^ , (1) 

/;<T 

where c^^, c;_cr are creation and annihilation operators, 
respectively, for electrons with spin a =t, i on site I with 
three-dimensional coordinate fi . The matrix elements ti 
are the electron transfer amplitudes between neighboring 
sites. Following Ref. 12 we set 

= 2.4494 eV , id = 2. 7939 eV , it = 3.4346 eV 

(2) 

for the single, double, and triple bonds, respectively. We 
consider the half-filled band exclusively, i.e., the number 
of TT electrons equals the number of lattice sites N . 
The electrical current operator is given by 

J = -iea ^ ti {c+_^^^„ci,a - c+^h+i.a) , (3) 

/,m;cr 

where a is an average bond length, i.e., we ignore geom- 
etry effects in J due to the difference in bond lengths 
Finally, we define the operator for the dipole moment, 

d-^lri-riKnj-l) . (4) 
I 
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Here ni = fii^i + hi^i counts the number of electrons on 
site I, and g- — c^^ci^cr is the local density operator at 
site I for spin a. Recall that we treat the PDA chain 
as perfectly straight, so that |r; — ri| is the appropriate 
sum over the bond distances rg, rj, and for the single, 
double, and triplet bonds between the sites / and m. 



2. Coulomb interaction 

The electrons interact electrostatically via the Cou- 
lomb interaction (Pariser-Parr-Pople modeli^) 



^ !/(/ - m) [{ni - 1) (n„, - 1)] . (5) 



The PDAs are insulators. Therefore, the Coulomb in- 
teraction is not dynamically screened at the energy scale 
of a few electron volts, and the screening is taken into 
account reasonably well by a static dielectric screening 
with dielectric constant Cd = 2.3 for PDAs. 

For the description of electrons and holes in quantum 
wires and other quasi one-dimensional structures, various 
effective potentials have been used in the literatureJ^ii^ 
For polymers, the general Pariser-Parr-Pople potential^*' 
is often approximated by the semi-empirical Ohno poten- 
tiaU^HOl 



V 



Ohno 



{I - m) 



V 



^l-f/3(|r1-r;„|/A)2 



(6) 



with U = V. At large distances, the Coulomb interaction 
must be recovered. Therefore, we require 



Vil- m) 



for |n-r„|>-^, (7) 



which implies y/J] = F/(14.397eV), where we have used 
that e^/(2aB) = 13.605 eV is the Rydberg energy and 
Ob — 0.5291 A is the Bohr radius. The remaining free 
parameter V describes the modification of the Coulomb 
potential at short distances due to the confinement of 
the electrons to the chain. Below, we derive the Ohno 
potential and justify it for intermediate to large length 
scales. For short distances, however, the Hubbard inter- 
action must be kept explicitly, leading to an additional 
parameter. 



3. Effective Coulomb potentials 

Our derivation of the one-dimensional effective po- 
tentials for the various cases closely follows Refs. [la 
and [3 In order to set up the single-particle basis in 
which the one-dimensional Hamiltonian ([5|) is formu- 
lated, we solve the single-particle Schrodinger equation 



for electrons whose motion is restricted to the 5;-direction 
due to a confining potential. We set Wconf (a^, y, z) = 
W2{x, y)Wi{z), so that the single-particle wave functions 
factorize: 'if{x,y,z) = ^{x,y)(j){z). The potential Wi{z) 
incorporates the (small) effects of the various ligands and 
permits a discrimination of the poly-diacetylenes. For the 
purpose of this work, we set it to a constant which we 
absorb into W2{x,y). 

The confining potential perpendicular to the chain di- 
rection is assumed to be very strong so that the electron 
wave function in the direction perpendicular to the chain 
is given by the lowest-energy state ^o{x, y)- The effective 
Coulomb potential between two charges at distance \z\is 
then given bj^i^^ 

V^^iz) . [,x,y,x',y'^MMmj:^ . 

J {x - x'Y + [y - y'Y + z"^ 

(8) 

This expression can be simplified further for a parabolic 
confining potential;^ 

1 9.9 9s '2h^{x'^ + y^) 
W2ix, y) = -mu:l,^,{x^ + y") = ' , (9) 

where the parameter Wconf — '2,Ti/{mE?) characterizes the 
strength of the confining potential. The ground-state 
wave function of the harmonic oscillators in the x and y 
directions is given by 



6(2;, y) 



2 1 



y 



i?2 



(10) 



where |'Co(2;,?/)P is a Gaussian with standard deviation 
Ax = Ay = R/2. This means that we find the elec- 
trons in the region (|a;| < i?, \y\ < R) with a probability 
of more than 90 percent. In poly-diacetylene single crys- 
tals, the distance between chains is typically d ~ 5 A so 
that the condition of weak overlap between the chains, 
i? ^ d, is fulfilled for R — 3.6 A. This view is sup- 
ported by the fact that the optical excitations for poly- 
diacetylene single crystals do not differ much from those 
for single chains diluted in their single-crystal monomer 
matrix. 2"'^ The excitation energy to the next confinement 
level is ?iwconf = 15.235 eV/(i?/A)2. For R = 3.6 A, 
the excitation energy to the next confinement level is 
hujconf — 1.2 eV so that higher confinement levels are not 
thermally populated at room temperature. 

When we insert (fTU]) into ([5]), we can carry out the 
Gaussian integrals in the coordinates X — {x + x')/2 and 
Y ~ {y + y')/2 and are left with a double integral over 
x^ = X — x' and y^ = y — y' ■ In polar coordinates, the 
resulting angular integral becomes trivial and we find^^ 



Jo 

2 

■v/tt exp 



dr(2r; 



exp {~{r/Rf 
\/ r"^ -\- z'^ 



e 

'r 



{z/RY [l-erf(|z|/i?)] ,(11) 
where erf (a;) is the error function. 
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We compare the long-distance limit of ([TTj) with that 
of the Ohno potential We demand the coefficients 
of the order l/\z\ and l/\z\^ agree and, with the help of 
Eq. (7.1.23) in Ref. 22, find that 



14.397 eV 

(R/ao) 



(12) 



or R/ao = with oq = 1 A the unit of length. In 

terms of the confinement parameter R, the screened po- 
tential and the Ohno potential can be cast into the form 



V 



Ohno 



(^) = 



V 



(13) 



y-'-(^) = V^V^exp[(z/i?)'] [l-erf(|z|/i?)],(14) 

whose large-distance expansions differ only to order 
{R/\z\)^. A comparison of the screened potential and 
the Ohno potential for all distances is shown in Fig. [51 
The agreement is to better than 10% for all \z\ > R. 
Even at \z\ = R/2 « 1.8 A, the discrepancy is only 25%. 
Therefore, it is justified to replace the screened poten- 
tial by the Ohno potential for intermediate distances. At 
short distances, the differences between the Ohno po- 
tential and the screened potential are substantial. The 
Ohno potential at z = is V^°^"°(0) = V, in contrast to 
V'^'iz = 0) = with ~ 1-77. 



1.5 
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B. Hamilton operator 

1. Hubbard-Ohno potential and screened potential 

In this work, we study the correlated motion of elec- 
trons along a chain, which we model using the Hamilto- 
nian 



H = T + V . 



(15) 



with T given by Eq. ^ and the potential V by Eq. ([S]). 
We use the Hubbard-Ohno potential 



v/l + (|n -7>„|/i?Ho)2 



(16) 



for / 7^ m and J7ho = V^Vho- Here the effec- 
tive Coulomb parameter Vho is linked to the screening 
length i?HO via Eq. Vho = 14.397 eV/(i?Ho/ao)- 

For example, we find Rao = 3.6 A for Vho = 4.0 eV. 
The Hubbard-Ohno potential was used earlier by Chan- 
dross, Mazumdar et al^ to explain the main absorption 
features of poly(para-phenylene vinylene) (PPV) with 
the parameter set Vcm = 4.0 eV and Ucm — 8.0 eV « 
\/t^Vcm- Our derivation in Sect. lII A SM ustifies their mod- 
ification of the Ohno potential, and their results for PPV 
support our choice for Vho = 4.0 eV. 

For comparison, we also give results for the screened 
potential. 



V^"{1 - m) 



VsciVttcxp 
1 - erf 



In - rrn\ 

RsCY 



In 



Rs, 



(17) 



which implies Usa- = -v/Tr^cr- Again, we have Vsa = 
14.397 eV/(i?scr/A). For example, we find i?scr = 4.1 A 
for Vscr — 3.5 eV. 



|z|/R 

FIG. 2: (Color online) Ohno potential V°''"°(z)/1/ from 
Eq. dlSj (full line) and screened potential {z)/V from 
Eq. (|14|l (dashed line) as a function of \z\/R. 



The above derivation applies to a straight geometry, 
not to the zigzag geometry of the Lewis structure in poly- 
diacetylenes. We account for the corresponding reduction 
of the chain length by the approximation |z| ~ |n — nnl 
as in the Ohno potential ^ , thereby ignoring the minor 
changes due to the non-orthogonality of the chain axis 
and the x-axis. 



S. Particle-hole symmetry 

The Hamiltonian ((T5)l and the current operator ([3]) are 
invariant under the particle-hole transformation c/,(t i— > 
(— l)'c^^. At half band-filHng, the ground state |$o) is 
also invariant under this transformation. Therefore, the 
expectation value of the dipole operator ([4]) vanishes in 

the ground state, do = ('i'oMl^o) = = 0. Likewise, 
the expectation values of the current operator vanishes 
in the ground state, (<i>o| Jj^o) — 0. 

We expect the same relation for excitons (bound 
particle- hole excitations of the ground state at half band- 
filling). Therefore, in the presence of a weak electri- 
cal field, such states show a quadratic Stark effect, see 
Sect. iHIBl 
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III. METHOD 



B. Polarizability and exciton wave function 



A. Single-particle gap and exciton binding energy 

Our model description contains a single parameter that 
can be taken to be the screening length R or the Coulomb 
parameter V . We use this parameter to adjust the typical 
band gap in poly-diacetylenes. The band gap or single- 
particle gap is defined by the difference in chemical po- 
tentials for a system with and A'c — 1 particles, 



= [Eo{N, + 1) - Eo{N,)] - [So(7Ve 



(18) 
1)] , 



which, due to particle- hole symmetry at half band-filling, 
reduces to 



A - 2 [Ea{N, ^N + 1)- EoiN, = N)] , 



(19) 



where Eo{Nc) is the energy of the iVo-particle ground 
state |$o)- Optical excitations to above the single- 
particle gap are extended and thus can transport cur- 
rent through the system. Experimentally, they can be 
monitored by the onset of the Franz-Keldysh oscillations 
in the electro- absorption signal. Typical values for poly- 
diacetylene single crystals are A = 2.3 eV for DCHD and 
A = 2.5 eV for PTS and PFBS.^a 

In poly-diacetylenes, the singlet exciton and its vi- 
bronic replicas carry most of the oscillator strength of 
the optical excitations. The quadratic Stark effect in the 
electro-absorption proves that they are bound states of 
electron-hole excitations The exciton energy thus de- 
fines the optical gap. 



^opt 



E,^iN, = N) - EoiN, ^ N) , 



(20) 



where Ec^{Nc — N) is the energy of the first excited state 
of the half-filled system |$cx), which has a finite overlap 
with an optical excitation of the ground state, J|<i>o)- 
The binding energy of the excitons is then obtained as 



A 



bind 



A, 



opt 



(21) 



For DCHD and PTS/PFBS PDA single crystals, the cor- 
responding binding energies are Abind = 0.48 eV and 
Abind = 0.51 eV, respectively!^ 

In our numerical investigation, we calculate the ground 
state |<I>o) and excited states |3>s), s > 1. An optical 
excitation with the current operator J, Eq. ([3]), has the 
oscillator strength 



< 1 



s>l 



Ws 



(22) 



Amongst the excited states l^s), we identify the exci- 
ton state l^'ex) as the energetically lowest- lying excita- 
tion that carries significant optical weight, > 0.1. 
The same description can be obtained using the dipole 
operator given by Eq. 



Since it is a bound state, the exciton displays a 
quadratic Stark effect, i.e., the redshift of the resonance 
level with an external static electrical field of strength F 
is given by 



SA. 



1 



opt 



-pF^ 



(23) 



where p is the polarizability. Note that the experiment 
measures the Stark shift both of the ground state |<I>o) 



mine Eo{F) and E^ 
Hamiltonian 



(F) for various fields F from the 



H{F) ^ H ~ eFd 



(24) 



with the dipole operator d, see Eq. (jH). Note that we 
determine the polarizability as measured in experiment, 
i.e., we need not resort to further theoretical considera- 
tions here. 

In order to extract the 'exciton radius' from the po- 
larizability, one can start from the Frenkel picture^^ or 
from the Wannier picture. As demonstrated in Ref. [25l . 
the probability distribution Pcx{l, m) provides a very de- 
tailed picture of the spatial character of the exciton in 
a many-particle approach. It describes the particle-hole 
content of the exciton wave function with respect to the 
ground state, i.e., it gives the probability that |$cx) is an 
electron- hole excitation of the ground state |$o) at sites I 
and m, respectively. Explicitly, 



Pex(^™) 



and 



(25) 
(26) 



We denote the probability density to find an electron- hole 
pair at a separation roh by Pcx(''oh), 

i^cx(0-I]^ex(;,m),5(r-|r5-r,„|) • (27) 

l,m 

The average electron-hole distance is then obtained 
from 



Teh = [r)cK = 



J drrPex(r) = ^Pc= 



{I, m)\fi - f„ 



In Ref. |20|, a simple two-level model was considered in 
which the exciton couples to a (representative) contin- 
uum state. As we shall see in Sect. IIVDI the two corre- 
sponding electron-hole distances compare well with each 
other. 
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C. Numerical procedure 

We present results both for the Hubbard-Ohno poten- 
tial, Eq. PB)) . and for the screened potential, Eq. ([TT]). 
In the presence of long-ranged Coulomb interactions, a 
high numerical accuracy is of crucial importance. There- 
fore, we devote this subsection to the problem of how we 
determine and control the accuracy of our calculations. 

In this work, we have performed the numerical cal- 
culations on finite chains with open boundary condition 
(OBC) using the non-local version^^'^^-^* of the DMRG 
techniquci^ The number of block states has been selected 
according to the dynamic block-state selection (DBSS) 
approach,22i^ 



1. Ground state and single-particle gaps 

In order to calculate the band and spin gaps, we have 
determined the lowest-lying eigenstates of various spin 
and charge sectors from independent DMRG runs. In this 
case, the DBSS approach-^^-'^" permits a rigorous control 
of the numerical accuracy because we can fix the thresh- 
old value of the quantum information loss x- Here we 
take X — 10~^. As another check, we have used the en- 
tropy sum rule for finite chain lengths for each DMRG 
sweep, i.e., we have verified that the sum rule has been 
satisfied after the third sweep. During our calculations, 
the maximum number of block states was varied in the 
range 256 < Mmax < 400 for OBC due to the large spin 
and band gaps. For M^ax = 400, the maximal chain 
length is A^max — 150. For these parameters, the individ- 
ual states can be treated reliably with an accuracy given 
by x= 10-5. 

In the presence of long-ranged interactions it is crucial 
to use a large Mmin in order to provide a good environ- 
ment block, i.e., to maximize the KuUback-Leibler en- 
tropy:— Here we take M^i^ — 128. Therefore, we have 
kept a number of block states that have small weight 
during the system build-up in the infinite-lattice step. 
During the sweeping iterations of the finite-lattice part 
of the DMRG, they gain more weight,— so that they sub- 
sequently become important for an accurate description 
of the ground state and the excited states. 



2. Optical excitations 

In order to calculate the optical gap, we have simulta- 
neously calculated low-lying eigenstates of the half- 
filled charge sector. In this case, it is necessary to spec- 
ify how the block entropy is calculated from the target 
states. In the presence of several target states, it is possi- 
ble to derive an upper and a lower bound for the mutual 
information between the system block and the environ- 
ment block. Therefore, an upper bound^ and a lower 
bound^ can be derived for the accessible information,^" 
but an exact expression is not available. 



In our work, we have defined the reduced density ma- 
trix before truncation as p = J^aPaPa, where the Pa are 
the reduced density matrices for the individual target 
states, and we have used fixed weights Pa > 0, J2a Pa — ^■ 
We have tried various sets of values for pa in order to sta- 
bilize the calculations and to improve the accuracy. For 
most of the results presented here, we have found that 
the statistically independent choice, Pa = l/Ag, provides 
satisfactory results. We have set the minimum number of 
block states to M^i^ — 400, and the maximal number of 
block states used in our calculations is Afmax = 800. We 
note in passing that the Davidson diagonalization rou- 
tine gives stable results only for < 2. For all cases 
of interest, A^s > 2, we have used the Lanczos method in 
order to obtain stable results. 

In order to identify the exciton state, we begin our cal- 
culations with 10 < A's < 15 target states with a reduced 
demand in accuracy. Once we have found the dominant 
optical excitation from the oscillator strengths p2)) . we 
repeat our calculations using the smaller number of tar- 
get states actually required, typically A^s = 5. We inde- 
pendently determine the exciton state from the optical 
weights based on the current operator, Eq. and the 
electrical dipole operator, Eq. (jl]). 



3. Chain topology 

In the standard DMRG procedure for OBC, it is more 
efficient to treat models which possess reflection symme- 
try; computational costs can be reduced significantly by 
applying the symmetry. In turn, the results are usually 
more accurate for the same parameter set and computer 
resources as compared to a non-reflection-symmetric con- 
figuration. In our study, a reflection-symmetric conflgu- 
ration can be realized by an appropriate choice of the 
bond sequence along the chain, e.g., [T(SDST)m] with 
N = 4m + 2 carbon atoms or [(SDST)™_iSDS] with 
A^ = 4m carbon atoms. This approach, however, has 
a drawback. We have found that, for the reflection- 
symmetric conflgurations, end excitations that are sim- 
ilar to the end spins for the S = 1 Heisenberg chain 
with OBC appear.— Due to these extra degrees of free- 
dom, the ground state becomes four-fold degenerate in 
the thermodynamic limit so that numerical calculations 
become less stable, especially when several target states 
are used to calculate the optical gap. In order to remove 
such end excitations, we could have modified the first and 
last spins or their couplings, as was done in Ref. fs^, or 
could have used the refiection symmetry as a conserved 
quantum number. In this work, we avoid these com- 
plications by using a chain configuration that does not 
have inversion symmetry, i.e., we use the bond sequence 
[(TSDS)™_iTSD] with A^ = 4m carbon atoms. 

In general, we have calculated the low-lying energy 
spectrum for both the symmetric and the non-symmetric 
chain configurations in order to identify the excitations 
unambiguously. In addition, we have determined the op- 
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tical gap and the dipole matrix elements for the current 
and dipole operators for both types of chain configura- 
tions. In the remainder of this paper, we present our 
results for the configuration [(TSDS)„_iTSD]. 



This limits the accuracy of the results obtained from 
finite-size extrapolations. In such cases, we restrict our 
extrapolations to use data up to A'^^ax = 80. 



4- Finite-size scaling 



IV. RESULTS 



The PDAs are charge and spin insulators, i.e., the gaps 
for single-particle, optical, and magnetic excitations are 
finite. The materials are characterized by finite correla- 
tion lengths. Therefore, end effects decay exponentially, 
and local operators that are calculated in the middle of 
the chain display a regular behavior as a function of in- 
verse system size. Thus, various quantities that we cal- 
culate for finite chain lengths N can be extrapolated re- 
liably to the thermodynamic limit, N ^ oo, by using a 
second-order polynomial fit. 
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FIG. 3: (Color online) Band gap as a function of inverse sys- 
tem size for various values of V for (a) the Hubbard- 
Ohno potential (|13|) and (b) the screened potential (1141) . The 
lines are quadratic fits. 



When we target several eigenstates simultaneously, our 
calculations on longer chains give less reliable results. 



A. Single-particle gap, optical gap, and exciton 
binding energy 

In Fig. [3] we show the single-particle gap as a func- 
tion of inverse system size 1/A^ for the Hubbard-Ohno 
potential p6|) and the screened potential (flT)) . The lines 
are quadratic fits in the inverse system size. The finite- 
size corrections to the result in the thermodynamic limit, 
TV ^ oo, are less than 0.05 eV for N > 100. 

As expected and as seen in Fig. [31 the single-particle 
gap increases as a function of the Coulomb parameter V. 
For the chosen band-structure parameters ([2]), the bare 
band gap is A{V = 0) = 1.25 eV, which is only half as 
large as the observed single-particle gap in PDAs. The 
Coulomb interaction accounts for the other half of the 
single-particle gap, i.e., exchange and correlations play 
an important role in this class of materials. In order 
to fit the experimentally observed gap, Acxp — 2.4 eV, 
we choose Vho = 4.0 eV for the Hubbard-Ohno potential 
and Vscr = 3.5 eV for the screened potential. These values 
correspond to a screening length of Ruo = 3.6 A and 
i?scr = 4.1 A, respectively. 
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FIG. 4: (Color online) Band gap and optical gap as a function 
of inverse system size 1/N for Vho ~ 4.0 eV for the Hubbard- 
Ohno potential (|16|) and for V^'^'^ = 3.5 eV for the screened 
potential (|17p . The lines are fits to a quadratic polynomial in 
1/iV. 

Both potentials display bound exciton states below the 
single-particle gap. In Fig. |3]we show the single-particle 
gap A (TV) and the optical gap Aopt(TV) as a function of 
inverse system size 1/N for both potentials. Figure 0] 
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shows that both effective potentials reproduce the exci- 
ton energy, Aopt « 1.9 eV, for PDA-DCHD. Correspond- 
ingly, we can reproduce the experimentally observed exci- 
ton binding energy, Abind = 0.5 eV, with both potentials. 



to detect the second exciton because it is hidden by the 
intense phonon replicas of the primary exciton. 



C. Exciton wave function and exciton radius 



B. Oscillator strengths, dark states, and second 
exciton state 

We show the distribution of oscillator strengths for the 
Hubbard-Ohno potential in Fig. [5l The screened poten- 
tial with Vscr = 3.5 eV leads to a qualitatively similar 
distribution in that the majority of the weight lies in 
the primary exciton. For = 80 sites, the oscillator 
strength for the primary exciton is Wex.HO = 0-59 for the 
Hubbard-Ohno potential for Vao = 4.0 eV. The exci- 
tons carry about 60% of the total spectral weight. This 
is in good agreement with experiment where the spectral 
weight of about = 1.2 of ntot = 2 7r-electrons is found 
below the single-particle gap; see Fig. 5 of Ref. [l^- 
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FIG. 5: Oscillator strengths Wb, Eq. (|22|l . as a function of 
energy for the Hubbard-Ohno potential (Vho ~ 4.0 eV) for 
N = 80. The weight of the primary exciton is ws — 0.59; 
the total weight of the first nine optically excited states is 
EL». = 0.60. 



We note that there are two optically dark states be- 
low the primary exciton, with energy differences A^i = 
0.24 eV and Aa2 — 0.05 eV. Both states are spin triplets. 
Experimentally, such dark states have been located at 
Aji = 0.4 eV."^— Therefore, we obtain a qualitatively and 
even semi-quantitatively correct ordering of the excited 
states. 

The primary exciton at the optical excitation energy 
Aopt(iV = 80) = 1.97 eV carries 99% of the excitonic 
weight. The second exciton around A'^p^{N = 80) = 
2.34 eV carries only 1% of the excitonic weight. Thus, our 
calculations indicate that two excitons should be visible 
in the PDA chains, whereby the second exciton has a 
binding energy of AJ^j^^^ = 0.1 eV and is lower in intensity 
by two orders of magnitude. Experimentally, it is difhcult 



In the following, we concentrate on the primary ex- 
citon. In Fig. [5] we show the probability distribution 
Pcxiljfn), Eq. ((26|) . i.e., the 'exciton wave function' for 
the Hubbard-Ohno potential for Vho = 4.0 eV. The 
probability distribution is similar for the screened po- 
tential for Vscr = 3.5 eV. 
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FIG. 6: Probability distribution Pcx(^,'ti) H26|) to find a hole 
at site I and an electron at site m in the exciton state. This 
'exciton wave function' is shown for the Hubbard-Ohno po- 
tential with Vho = 4.0 eV on a chain with A = 80 sites. 



The probability distribution reflects the structure of 
the unit cell with four carbon atoms. We expect that the 
exciton wave function factorizes, Pex(^TO) ~ ^'cm[(^ + 
m)/2]ip{\l — m\), where 5'cm[(^ + w)/2] describes the mo- 
tion of the center-of-mass and (p{\l — m\) describes its 
internal structure. The center-of-mass wave function fol- 
lows that of a particle in a box. In an infinite system, it 
corresponds to a state with zero total momentum because 
the light field adds only a negligible momentum to the 
ground state. In fact, for fixed |Z— m|, we observe nodes at 
the boundaries (/, m — > 0, A^) and a maximum in the mid- 
dle of the chain for m « A^/2. For fixed center-of-mass 
coordinate I -\- m, the probability distribution Pey^{l,m) 
reveals the internal structure of the exciton, (p{\l — m|). 
Cross sections of Pex(^ m) along the lines I + m — const, 
show that electron and hole are bound to each other, i.e., 
Lp{\l — m\) is vanishingly small for |Z — to| > r,„ax. 

The exciton wave function Po-x{l,m) shows a promi- 
nent odd-even effect as a function of the particle-hole 
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separation \l — m|. This is a consequence of the invari- 
ance of the Hamiltonian and of the current operator un- 
der a particle-hole transformation, see Sect. IIIB 21 At 
half band- filling, the ground state |<i>o) is invariant un- 
der this transformation. If the same applies to an ex- 
cited state an inversion symmetric system obeys 
{^s\c+Jm.a\'^o) = (-l)'+'"+i($.|c+^c„,,|$o). There- 
fore, the overlap vanishes for even \l — m\. Since the 

exciton obeys |$ex) = y/w^J\^o) + I'i'rcst), and the sys- 
tem is approximately inversion symmetric, it is only the 
(small) contribution |$rcst) which contributes to the ex- 
citon wave function for even I / — ml. 




FIG. 7: Radial distribution function Pcx(r) (l27ll for the 
Hubbard-Ohno potential ( Vko = 4.0 eV) for a chain of 
iV = 80 sites. 



Further insight into the structure of the exciton wave 
function is gained from the radial distribution function 
Pcx{r), Eq. (|27|) . for the electron-hole distance, which is 
shown in Fig. [7] for a chain of = 80 sites. As expected 
for a bound electron-hole pair, the distribution decays 
rapidly as a function of the electron-hole distance r. The 
peak at a distance r ^ 4 reflects the fact that there are 
four carbon atoms in the unit cell. The oscillations in the 
radial probability are the result of the odd-even effect ob- 
served in the probability distribution function Pcx(', fn), 
Eq. (pS)) . The overall behavior of the radial distribution 
resembles the results obtained from the GW^+BSE ap- 
proach to polymers j^iS, 

Finally, we show the exciton radius, Eq. as a func- 
tion of system size for the Hubbard-Ohno potential for 
Vho — 4.0 eV and the screened potential for K,cr = 3.5 eV 
in Fig.[8l As expected for a bound exciton, it does not in- 
crease much with system size for N ^ AO. We find an ex- 
trapolated exciton radius rf^ = 9.67 A and r^" = 8.54 A 
for the Hubbard-Ohno and screened potentials, respec- 
tively. 

For completeness, in Fig. [SI we show the wave function 
of the second, weak singlet exciton. The center-of-mass 
coordinate again describes a particle in the box, whereas 
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FIG. 8: Exciton radius as a function of system size for the 
Hubbard-Ohno potential (Vho ~ 4.0 eV) and the screened 
potential [Vscv — 3.5 eV). 



the internal structure for the relative motion of electron 
and hole displays a node as a function of \l — m\ for 
I + m =const. The size of the second exciton is about a 
factor of two larger than the size of the primary exciton. 
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FIG. 9: Probability distribution Pc^{l,m) (|26|) for the second 
exciton for the Hubbard-Ohno potential with Vho = 4.0 eV 
on a chain with TV = 80 sites. 



D. Polar izability 

The polarizability follows from the Stark shift of the 
exciton in the presence of an external electric field. When 
we measure the strength of the electrical field F in terms 
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of the energy unit / = Feao, with «□ = 1 A, we can write 
the polarizabihty in the form 



^ = 28.694-A_^ , 
ao-^ (eV) 1 



(29) 



where A describes the excitonic Stark shift, <5Aopt 
-A/2, see Eq. 
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FIG. 10: (Color online) Stark shift of the exciton binding 
energy as a function of the rescaled electrical field / = Feao 
(ao = lA) for various system sizes for the Hubbard-Ohno 
potential (Vho = 4.0 eV). 



Fig. [To] shows the redshift of the binding energy for 
the primary exciton due to the electric field for various 
system sizes. For systems N > 40, the curvature A does 
not change significantly. This reflects the fact that the 
exciton wave function does not depend on the system size 
for > 40. Taking a rough typical value from the fits 
in Fig. [TUl we estimate A w 3.0 • 10^ (eV)"^ so that we 

find p = 8.6 • 10"^ A'^ for the polarizabihty. This compares 
favorably with the experimental value for PDA-DCHD, 
PDCHD = 8.2 • 103A^ or PDA-PTS, ppts = 7.2 • lO^A^^ 
In the experimental workj^S a semi-empirical model 
was used to extract the exciton radius rexc from the po- 
larizabihty, p = a(ercxc)^/Aopt, where a is a factor of 
the order of unity. For a — 1, the experimental value for 
p leads to Tcxc ~ 12 A, in good agreement with our value 
for the average electron-hole separation, r^h = 9.7 A. 



E. Triplet exciton 

Finally, we summarize our results for the triplet sec- 
tor, i.e., excitations with total spin S — I. Note that it 
is difficult to access this spin sector experimentally. In 
Fig. [Til we show the differences in the ground state en- 
ergies in the spin singlet sector {S = 0) and the spin 
triplet sector (5' = 1), together with the optical gap in 
the triplet sector as a function of inverse system size for 



the Hubbard-Ohno potential (Vho = 4.0 eV). As can be 
seen, one finds two finite gaps with different sizes which 
we will interpret below. The results for the screened po- 
tential (Vscr = 3.5 eV) are very similar. 
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FIG. 11: (a) Energy difference between the ground states of 
the two spin sectors and (b) the energy of the triplet exciton 
as a a function of inverse system size 1/A'' for the Hubbard- 
Ohno potential (Vho = 4.0 eV). 



In Fig. [13 we depict the resulting energy level scheme 
for the Hubbard-Ohno potential for V = 4.0 eV. The 
energies are extrapolated values in the thermodynamic 
limit. The ground state of the triplet sector is Eq^^ — 
E^=° = 1.7eV higher in energy than the ground state for 
the singlet sector, see Fig. [TTk . Therefore, this state is 
frequently called the 'triplet exciton'. This state is 0.2 eV 
below the singlet exciton. Experimentally, however, this 
state has been detected at 0.9 eV below the singlet exci- 
ton. Therefore, our description underestimates this en- 
ergy, i.e., the triplet ground state should be much lower 
in energy. This discrepancy could be the consequence 
of the large polaronic effects in the triplet sector, as has 
been suggested in Ref. [H. The dominant optical excita- 
tion in the triplet sector lies only 0.6 eV above the triplet 
ground state, see Fig. [TTb. but still 0.4 eV above the sin- 
glet exciton and is close to the single-particle gap, i.e.. 
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FIG. 12: Energy level scheme for the spin singlet and 
spin triplet sector for the Hubbard-Ohno potential (Vho = 
4.0 eV). 



only 0.1 eV below the particle-hole continuum. 

The wave function of the triplet excitations is shown 
in Fig. [T31 Similarly to the singlet exciton discussed 
in Sect. IIV"Cl the probability distribution approximately 
factorizes into a center-of-mass wave function and a wave 
function for the relative coordinate. The center-of-mass 
wave function are similar in all cases. The wave func- 
tion for the relative motion of electron and hole shows 
two nodes for the optical excitation of the triplet ground 
state, whereas it has one node in the case of the optical 
excitation of the singlet ground state. 



CONCLUSIONS 
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FIG. 13: Probability distributions Pex(Z,m), Eq. ([261), for (a) 
the triplet ground state and (b) its optical excitation for the 
Hubbard-Ohno potential with Vho — 4.0 eV on a chain with 
iV = 80 sites. 



In this paper, we have studied the role of electron- 
electron interactions in poly-diacetylenes. Since the bare 
band gap is only half as large as the observed single- 
particle gap and since the binding energy of the singlet 
exciton of 0.5 eV is 20% of the single-particle gap, ex- 
change and correlation must play an important role in 
this class of materials. Our density-matrix renormaliza- 
tion group method permits the numerically exact treat- 
ment of an appropriate model Hamiltonian with long- 
range Coulomb interactions for a large number of elec- 
trons. We have used the experimentally observed single- 
particle gap, Aband — 2.4 eV, to fix the strength of the 
Coulomb interaction at short distances. We have deter- 
mined the lowest-lying optical excitation and have re- 
produced the observed exciton binding energy in poly- 
diacetylenes, Abind — 0.5 eV. 

The key difference between our work and previous nu- 
merical DMRG studies^ii^ is the parameterization of 
the Coulomb interaction. We argue that the Ohno pa- 
rameterization of the Pariser-Parr-Pople interactionl*^ is 
not appropriate for very short distances because the lo- 



cal (Hubbard) interaction is too small. We propose to 
use the Hubbard-Ohno potential fTB]) or the full expres- 
sion P?)) of the Coulomb interaction for effectively one- 
dimensional structures. Our findings support earlier the- 
oretical studies of poly(para-phcnylcne vinylene)i^ 

Our results indicate that the PDA chain has two opti- 
cally dark states below the exciton, in qualitative agree- 
ment with experiment. Moreover, the screening poten- 
tial supports a second bound exciton with binding energy 
^bind — 0.1 eV whose intensity is two orders of magni- 
tude smaller than that of the primary exciton. 

The exciton wave function approximately factorizes 
into two terms: the center-of-mass wave function and 
the relative wave function. The former describes the ex- 
citonic 'particle-in-a-box' state; the latter describes the 
pair state of an electron and a hole whose separation 
(exciton radius) rapidly converges to a finite value with 
increasing system size. In order to investigate the behav- 
ior of the exciton in the presence of an external electric 
field, we have studied the polarizability due to the Stark 
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shift of the exciton. As expected for a bound electron- 
hole pair, the exciton displays a quadratic Stark redshift 
in energy as a function of the field strength. Our calcu- 
lated polarizability reproduces the experimental results 
for PDA-DCHD and PDA-PTS chains. 

Finally, we have studied the triplet sector. The energy 
of the triplet ground state found in our calculation is too 
high, i.e., the binding energy of the lowest-lying triplet 
excitation is too small. We attribute this discrepancy 
to strong polaronic effects in the triplet sector In our 
work, we have not considered the effects of lattice relax- 
ation, the electrostatic potential of the poly-diacetylene 
side-groups, and geometry effects. The inclusion of these 
effects is required for a more detailed description of indi- 
vidual members of the poly-diacetylene family. 
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